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ABSTRACT 



Hydraulic ram concerns the dynamic loading and cata- 
strophic failure of liquid filled fuel tanks impacted by high 
speed projectiles. Hydraulic ram is divided into two phases; 
shock phase and cavity phase. The shock phase was studied. 
Theoretical predictions of shock radius versus time were com- 
pared with experimental results. The theory was found to be 
reasonably accurate for blunt shaped projectiles. 

Included in the theory were predictions of shock pressure 
and pressure profiles behind the shock. Theoretical pres- 
sures are presented graphically but were not validated by 
experiment. 

Velocity decay of the projectiles in the liquid was in- 
vestigated and correlated with theory to provide information 
for further studies of the cavity phase. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION 8 

II. ANALYTICAL STUDY 11 

A. SHOCK FRONT PRESSURES 12 

B. PRESSURE FIELD BEHIND THE SHOCK 13 

C. SHOCK MACH NUMBER 15 

III. EXPERIMENTAL PROCEDURE 17 

IV. RESULTS AND CONCLUSIONS 31 

A. SHOCK RADIUS VERSUS TIME - PRE-PUNCHED FRONT 

WALL 31 

B. SHOCK RADIUS VERSUS TIME - SOLID FRONT WALL 37 

C. PROJECTILE VELOCITY DECAY IN THE FLUID 40 

D. PRESSURE VERSUS RADIUS AND TIME 49 

V. RECOMMENDATIONS 59 

COMPUTER PROGRAM 61 

LIST OF REFERENCES r 67 

INITIAL DISTRIBUTION LIST 68 

FORM DD 1473 69 



3 



LIST OF FIGURES 



1. Ballistic Range Components 18 

2. Ballistic Range (Down-range View) 19 

3. Projectile Shapes 20 

4. Projectile Parameters ' 20 

5. Adjustable Time Delay Unit Circuitry 22 

6. Schematic of Shadowgraph Apparatus (Down-range 

View) 23 

7. Shadowbox, Tank and Stand Installation 25 

8. Front Wall Test Specimen, Pre-punched Hole 26 

9. Shadowgraphs After Projectile Impact Through a 

Pre-punched Front Wall 29 

10. Shock Radius vs. Time, Eq = 1600 in-lb 32 

11. Shock Radius vs. Time, E^ = 7493 in-lb 33 

12. Shock Radius vs. Time, E^ = 12323 in-lb 34 

13. Shock Radius vs. Impact Energy, E^ 35 

14. Solid Front Wall After Test 38 

15. Typical Damage to Solid Front Wall at Point of 

Impact 38 

16. Damage to Solid Front Wall (Mounted on Tank) 39 

17. Qualitative View of Projectile Energy and Velocity 

Decay in the Fluid 42 

18. Typical Deformation of Projectiles After Coming to 

Rest in the Fluid Filled Tank - Pre-punched Hole 
Results 4 4 



4 



19. Projectile Distance in Fluid vs. Time, Eq = 1600 

in-lb 45 

20. Projectile Distance in Fluid vs. Time, Eq = 7493 

in-lb 46 

21. Projectile Distance in Fluid vs. Time, E^ = 12323 

in-lb 47 

22. Pressure vs. Radius from Impact - Complete Profile, 

Eq = 1600 in-lb 50 

23. Pressure vs. Radius from Impact, E^ = 1600 in-lb — 51 

24. Pressure vs. Radius from Impact, Eq = 7493 in-lb — 52 

25. Pressure vs. Radius from Impact, Eq = 12323 in-lb - 53 

26. Pressure vs. Time for Three Radii, Eq = 1600 

in-lb 55 

27. Pressure vs. Time for Three Radii, E =7493 

o 

in-lb ■ 56 

28. Pressure vs. Time for Three Radii, E^ = 12323 

in-lb 57 



5 



LIST OF SYMBOLS 



A Area of deformed surface of the projectile 

B Constant in the equation of state; constant in velocity 
decay equation 

Speed of sound in undisturbed fluid 
Cjj Drag coefficient 

D Drag acting on projectile 

Eq Initial kinetic energy of the projectile; energy re- 
leased at a point source 

AEg Change in projectile energy due to shock formation 
H Constant defined by equation (2.14) 

Mg Shock Mach number 

N Exponent defined by equation (2.14) 

P Pressure 

P Pressure in undisturbed fluid 

o 

Rg Shock radius 

Rg Shock velocity 

Rg Shock acceleration 

Rq Initial radius defined by equation (2.17) 

T Temperature 

V Projectile velocity 

Projectile impact velocity 

AV’ Change in projectile velocity due to shock formation 

W 

e Specific internal energy of the fluid 

Specific internal energy of undisturbed fluid 
e^^ Specific internal energy after the shock 



6 



f Non-dimensional pressure; equivalent flat plate area 
Non-dimensional pressure after the shock 
g Non-dimensional internal energy 

g^ Non-dimensional internal energy of undisturbed fluid 

gj^ Non-dimensional internal energy after the shock 
m Mass of projectile 

n Exponent defined in the equation of state 

q Exponent defined in density profile equation 

r Radius, radial coordinate 

s Specific entropy of fluid 

t Time 

u Particle velocity 

V Specific volume (1/p) 

X Projectile distance traveled in tank 
3 Shock decay coefficient 

5 Non-dimensional radius 

p Density of the fluid 

Pq Density of the undisturbed fluid 

P£ Density of the fluid 

4) Non-dimensional particle velocity 

(|)j Non-dimensional particle velocity after the shock 

Non-dimensional density 
4^^ Non-dimensional density after the shock 



7 



I . INTRODUCTION 



Non-nuclear survivability of aircraft flying in a hos- 
tile arena has been studied throughout the history of avia- 
tion. Recent experiences in Vietnam have demonstrated the 
vulnerability of sophisticated aircraft to both small arms 
ground fire and the more complex surface-to-air defenses. 
Except for the rare case of a direct impact by a surface- 
to-air missile, the damage has been incurred by relatively 
small projectiles. 

Tactical aircraft fuel tanks or cells have the largest 
surface area and volume of all vulnerable components. For 
this reason fuel cells have been most susceptible to damage. 
Projectile impact with fuel cells have caused fuel starva- 
tion, fire or explosion as well as extensive structural dam- 
age. The survivability of aircraft when a ballistic threat 
causes destruction of fuel filled tanks has been the main 
impetus behind this study. Astronomical costs of combat 
aircraft have made it cost effective to acquire more knowl- 
edge, of the destructive mechanism to fuel cells. Another 
application of this information could be in space flight 
where meteoroid collisions are possible. 

Projectiles which penetrate fluid filled cells cause 
damage many times more severe than that incurred by impact 
with an empty cell. With sufficient projectile energy there 
will be catastrophic failure of cell structural components. 
Hydraulic ram is the term used to describe such destructive 
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events. The hydraulic ram phenomenon is usually described 
by two phases: the shock phase and cavity phase. The effects 
of these phases are assumed to be sequential, occurring in 
the order listed above. Upon initial projectile impact, 
energy is transferred to the fluid, creating a strong hemi- 
spherical shock centered at the point of impact. The shock 
pressure loading may cause failure of the entrance wall. 
Pressures load the wall for a very short time because of 
geometrical attenuation. Essentially there is an impulsive 
load on the front wall during the shock phase. 

As the projectile travels through the fluid, a cavity is 
formed. This phase has sometimes been called the bubble 
phase because of the bubble-like appearance of the cavity. 
Others have called this phase the drag phase, using the la- 
bel to describe the mechanism which creates the cavity [Ref. 
11 . Oscillation of the cavity volume produces pressure 
pulses that may be sufficiently intense to cause catastroph- 
ic failure of cell components. The intensity of the result- 
ing pressure pulses is weaker than the shock phase pressure 
pulse but the duration is much longer [Ref. 2] . Ultimate 
destruction of tank components other than the entry wall has 
been generally attributed to the cavity phase. 

Until recently, experimental data from studies of hy- 
draulic ram have not been able to describe the basic phenom- 
enon. Entire aircraft or wing sections have been destroyed 
while studying the phenomenon. This method of experimenta- 
tion has become too costly; hence theory explaining the 



9 



phenomenon is necessary to defeat hydraulic ram. Much of 
the theory as it applies to hydraulic ram is still in its 
infancy. The complete mechanism is complex, requiring ac- 
curate and detailed analytical efforts. Correlation of 
theory with experiments must be accomplished simultaneously 
to assess theory validity. 

Research at the Naval Postgraduate School has been de- 
signed to isolate and observe the individual aspects of hy- 
draulic ram. The ballistic range has been used to study 
energy losses through aluminum plates, shock propagation and 
the cavity phase. Specific details were observed with the 
intent of applying the knowledge gained to the design of a 
survivable fuel cell. 

This study has investigated the shock phase using the 
theory of Yurchovitch and Lundstrom [Refs. 1 and 2] . A com- 
puter program was formulated from the mathematical model. 
Experimental data for various projectile energies were com- 
pared with the results of the theory. Yurchovitch' s formu- 
lation of shock radius as a function of time has agreed with 
the experimental data. These tests also have shown that 
damage to the tank front wall was not extensive enough for 
significant fuel leakage at the energy levels tested. 
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II. ANALYTICAL STUDY 



In order to predict the effects of the shock phase of 
hydraulic ram, it is necessary to study how the projectile 
impact energy affects the fluid that is disturbed. It was 
assumed that the impact energy, E^, was known and was suffi- 
cient to create a strong hemispherical shock. The perturba- 
tion in the fluid can then be described by assuming that the 
shock radius is proportional to time to the 0.8 power, and 
that there is a power law density profile behind the shock. 
This formulation was suggested by Yurchovitch [Ref. 2] . 
Yurchovitch found that for a given fluid equation of state 
the energy transferred to the fluid determines the resulting 
flow field. 

A. SHOCK FRONT PRESSURES 

The one-dimensional Rankine-Hugoniot equation for a 
shock moving through undisturbed fluid was used' to determine 
shock front conditions. To complete the solution an equation 
of state was needed. The Tait Equation for isothermal com- 
pressibility of liquids, 

- i(-^) = i , (2.1) 

V 3Pg n [P+B (s) ] 

where n is a constant for a specific liquid and B is a func- 
tion of entropy only, was integrated along lines of constant 
pressure and entropy. Integrating equation (2.1) yields 
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( 2 . 2 ) 



P 



B(S) [( 



V(T,0) 

V(T,P) 



n 

) - 1 ] 



Assuming the changes across the shock are adiabatic and the 
pressure in the undisturbed fluid is small compared to the 
shock front pressure, B is given as 

B = where = (^) (2.3) 

n ° s 

With the assumption that B and n are constant for a given 
fluid the equation of state is 

P = •^^[(— )^-l] (2.4) 

n Pq 

For water n, a constant, was found to be 7.0 [Ref. 3]. 

Non-dimensionalization of the variables was found to be 
convenient and the following scheme was adopted: 
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Ms = Rs't)/Co 
l[>(C,Mg) = 



<[>(?, Mg) 



u(r,t) 

ftg(t) 



f (?,Mg) 



p(r,t) 
Po^s (t) 



g(C,Mg) 



e (r, t) 
R|(t) 



(2.5) 
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With the assumption that Pq<<P and eQ<<e the Rankine-Hugoniot 
equations become 



1 = 

= <fi 



( 2 . 6 ) 




and the equation of state is 






(2.7) 



nM^ 



Solving the above equations for the parameters at the shock 
the following is obtained : 



By assuming values for the equation (2.8) can be solved 

for the shock Mach number. Mg. Subsequently, fj and ({> j ver- 
sus Mg can be found. Conditions at the shock front are now 
known while conditions behind the shock front are still un- 
known. 

B. PRESSURE FIELD BEHIND THE SHOCK 

The density profile behind the shock is assumed to have 
the form 



1 



( 2 . 8 ) 





(2.9) 



( 2 . 10 ) 






( 2 . 11 ) 
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To find q(M ) the total mass enclosed by the shock/ 

2d 

PoR|2TT/ii;5M5 , 

is equated to the original mass occupied by the enclosed 
volxome/ 

2/3ttR|Po 
with the result 

= 1/3 (2.12) 

Q 

Substituting the assumed density profile into (2.12) > q(Mg) 
is found to be 

q(Mg) = 3(t|)j-D (2.13) 

and the density function becomes 



It is now possible to solve for the particle velocity, <p 
by substituting equation (2.12) into the one dimensional con 
tinuity equation and solving for the constant with (()=(}) j , at 
5=1 (the shock front location) . The continuity equation is 
given by 



35 35 5 ®3M, 



where 3 = R_R_/R^ is the shock decay coefficient. With the 
above substitution and integration, the particle velocity 



profile is found to be 

<P = <i>iC 



By substituting the density and particle velocity pro- 



files into the momentum equation 




s 



and integrating, the pressure profile is found to be 




where the constant of integration was solved with f=f j at 



The shock- decay coefficient, 6, can be found using the 
assiamption that shock radius is given by 



Experimental work done by NASA [Ref. 4] has found that 



for water. Substituting equation (2.14) into the expression 
for 3/ the decay coefficient becomes 
3 = -0.25 

It is now possible to solve the pressure distribution 
function given the shock Mach number. Mg. 

C. SHOCK MACH NUMBER 

By considering the initial energy being released at the 
center of symmetry in the fluid it is possible to equate the 
initial energy to the integral energy density behind the 
shock. The energy liberated initially, Eq, is distributed 
such that 



R. 



s 



(2.14) 



N = 0.8 




(2.15) 
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Using the assumptions of a strong hemispherical shock as well 
as the density profile and particle velocity profile as de- 
rived previously, the integral equation (2.15) in non-dimen- 
sional form becomes 



Eo = (2.16) 

The values of p^,.C^ and n can be found for a given fluid. 
By choosing a Mach number the dependent variables (})j and 
are known. Equation (2.16) can then be integrated and 
2 Eq/R| versus shock Mach number can be plotted. Given the 
energy released, E^, it is then possible to find the shock 
Mach number as a function of shock radius from the plot of 
2 E_/r 1 versus M . 

There is however, an anomaly at a shock radius of zero. 
To avoid a finite amount of energy in zero volume or an in- 
finite energy density, reference 5 suggests the initial 
radius be given as 



= (- 



Po^o^'^ 



-) 



(2.17) 



Equation (2.17) was then used as the starting point for the 
solution of the theoretical pressures behind the shock. 
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Ill, EXPERIMENTAL PROCEDURE 



The basic components of the ballistic range used for this 
study are shown in figure 1 . A down-range view of the bal- 
listic range in figure 2 depicts the chronograph screen ar- 
rangement and test tank. The two spark sources above the 
tank were used as light .sources for the shadowgraphs. Two 
rifles were used, a .22 caliber and a Remington .222 caliber. 
Commercial Remington .22 caliber Long Rifle ammunition was 
used. The .222 ammunition was hand loaded for testing so 
that the velocities of the projectiles could be adjusted and 
would be known more precisely. Sierra bullets classified as 
spitzers were used. Figure 3 displays the corresponding 
shapes of the projectiles. The mass of the .222 caliber and 
.22 caliber projectiles are listed in figure 4 along with 
the parameters used to calculate the energy of each configu- 
ration. The area and drag coefficient parameters listed in 
figure 4 were experimental values used to calculate the ve- 
locity decay that was experienced during projectile passage 
through the test tank. 

AVTRON, No. A914T333, chronograph screens were used to 
provide start and stop pulses to Monsanto 101 B counters. 

As the projectile passed through the screen the pulse circuit 
was energized. The first screen started the first counter 
and the second screen stopped the first counter and simultan- 
eously initiated the second counter. The third screen 
stopped the second counter. The projectile impact velocity 

was computed using these times and distance measurements. 
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Figure I. Ballistic Range Components 
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Figure 2. Ballistic Range (Down- range View) 




Figure 3, Projectile Shapes (L.to R.; .22 
caliber, m=5.75 x 10“^ lb, .222 caliber, m= 6.43 
x|0"3|b, .222 caliber, m = 7.85xl0"3 lb) 



Caliber 


.22 Long Rifle 


.222 Remington 


.222 Remington 


Impact Velocity, Vo 
(ft /sec) 


1200 


2500 


2900 


Projectile Mass, m 
(Ibm ) 


5,75xI0‘3 


6.43x10*^ 


7.85x10"^ 


Impact Energy, Eq 
( in-lb) 


1600 


7493 


12323 


Experimental Drag 
Coefficient, Cq 


1.20 


0.45 


0.50 


Frontal Area of 
Projectile After 
Traversing Fluid, A 
(in2) 


0.122 


0.098 


0.165 


Equivalent Flat 
Plate Area, f(=CrvA) 

Ciu^) 


0.146 


0.044 


0.083 



Figure 4. Projectile Parameters 
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The events within the fuel cell were timed with two addi- 
_tional counters. A fourth chronograph screen provided a 
signal pulse which was distributed to the adjustable time 
delay unit (ATDU) and the two remaining counters to initiate 
them. The main element of the adjustable time delay unit 
circuitry was a Signetic, N74123, retriggerable, monostable 
multivibrator as shown in figure 5. Within the ATDU were 
two duplicate delay circuits to provide a 100 volt output 
pulse to each spark source. The time delay for both circuits 
in the unit was adjustable; however, the duplicate circuit 
delay was initiated by the first circuit's output. This re- 
sulted in the delay time for the second output pulse from 
the ATDU being the sum of the delay time of the first cir- 
cuit plus the delay of the second circuit. The delay time 

range for the first circuit was 500 to 1600 micro-seconds, 
while the second circuit provided an additional delay of 5 
to 15 micro-seconds. The 100 volt output pulses going to 
the spark sources were also used to stop the two counters. 

The first output pulse went to the spark source providing 

illumination for a shadowgraph in the horizontal plane, 
while the second pulse triggered the spark source for a 
shadowgraph in the vertical plane. Each counter therefore 
measured the time the shadowgraphs were taken. 

Two 16-inch collimating mirrors were used to pass light 
through the tank. The spark sources were located at the 
focal point of their respective mirrors. Figure 6 demon- 
strates the arrangement of the tank with respect to the 
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Vcc ( 5 Volts d.c.) 




Figure 5. Adjustable Time Delay Unit Circuitry 
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Figure 6. Schematic of Shadowgraph Apparatus 
(Down-range View) 
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spark sources, mirrors and photographic plates. The side 
and bottom walls of the tank were made of plexiglass through 
which very little light attenuation occurred. One side wall 
and the bottom wall had grid lines etched in the plexiglass 
for use in measuring the shock radius and cavity distance. 
High speed Polaroid film. Type 57, was used for the shadow- 
graphs. A ground, first surface mirror was used to deflect 
the light in the vertical plane 90 degrees. 

The tank or test cell had front wall inner dimensions of 
17.5 inches wide and 17.0 inches high with a depth of 17.0 
inches. Figure 7 contains a view of the tank. The main mem- 
ber of the front wall was one-inch thick aluminum and served 
as a mounting plate for the specimens used to simulate fuel 
cell wall material. To accommodate the test material the 
mounting plate had a circular cut-out 15.0 inches in diameter. 
The specimen was mounted behind the mounting plate with 16 
bolts arranged in a circular pattern. A front wall of 7075- 
T6 Aluminum, .090 inches thick, with a prepunched entry hole 
is shown in figure 8. This wall was used in the experiments 
with the exception of the few tests which used a solid front 
wall. The rear wall of the tank was identical to the front 
wall. 

To provide support for the front and rear walls aluminum 
angle strips were placed longitudinally at the four corners 
of the tank. The angle strips were bolted to the front and 
rear mounting plates. In addition to providing structural 
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Figure 7. Shadowbox, Tank and Stand Installation 
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Figure 8. Front Wall Test Specimen, 
Pre-punched Hole 
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rigidity, the angle strips retained the one-inch thick plexi- 
glass tank walls. The top of the tank was open. 

Room Temperature Vulcanizing (RTV) rubber sealant was 
used at all joints to make the tank water tight. At the two 
highest energy levels the sealant around the front wall had 
to be replaced several times. Apparently the shock pres- 
sures were sufficient to deteriorate the sealant bond to the 
plexiglass and aluminum. Tape was used to prevent water 
spillage through the prepunched hole and to avoid any signif- 
icant energy loss during penetration. 

During initial testing the tank was set on the top of 
the stand. It was discovered that although the tank was re- 
strained by the stand in longitudinal and lateral motion, 
angular tank movement was created by the moment arm and the 
force of impact at the center of the tank. This invalidated 
the shock radius data and the tank was subsequently bolted 
to the frame. It was also discovered that the tank and 
stand would move approximately eight inches during a test 
firing. To alleviate this problem the front two legs of the 
stand were bolted to the floor. Eventually these bolts 
worked loose and finally all four legs were bolted to the 
floor and 150 pounds of sandbags were placed on the stand. 

The experiment was run at three energy levels. One en- 
ergy level was tested before continuing with the next higher 
energy level. By varying the time delay with the ATDU, 
shadowgraphs of the shock radius versus time were taken. It 
was also possible to measure the projectile distance into 
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the tank versus time. Figure 9 shows a pair of typical 
shadowgraphs. The impact velocity of the projectile was 
computed from the times on the first two counters and the 
distance between chronograph screens. The shadowgraph times 
were measured by the third and fourth counters. 

Acquiring data for the .22 caliber was found to be very 
difficult because of the deviation in impact velocities from 
one test firing to another. A good estimation of velocity 
was required prior to a test to compute projectile time of 
flight from the fourth chronograph screen to the tank wall. 
That time plus the desired time for photographing the shock 
radius determined the delay time used to set the adjustable 
time delay circuit. Since the .222 caliber projectiles were 
hand loaded, the impact velocities were more predictable. 
However, in the .222 caliber tests the velocity had a suffi- 
cient deviation to make it difficult to shadowgraph a shock 
radius smaller than one inch. 

The data were reduced by first calculating the projectile 
impact velocity. The time of flight of the projectile from 
the fourth chronograph screen to the tank was then calcu- 
lated. This time was subtracted from the elapsed time of the 
shadowgraphs as measured by the third and fourth counters. 

The result was the time corresponding to the shock position 
from initial impact at the front wall. Shock radius, meas- 
ured on the shadowgraph, v/as then plotted at the correspond- 
ing time. Measurements of the shock radius were accurate to 
±0.02 inches. 
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Figure 9. Shadowgraphs After Projectile Impact Through 
a Pre-punched Front Wall 
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Projectile distance into the tank versus time data were 
also acc\amulated. The shadowgraph at 24 micro-seconds in 
figure 9 shows that the projectile was not distinctly visible. 
Work done by McMillen [Ref. 6] implies that the projectile 
nose is at the tip of the light cusp seen at the leading edge 
of the cavity shadow. Light refraction caused by density 
gradients in the flow field forms this nearly triangular re- 
gion called the cusp. 
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IV. RESULTS AND CONCLUSIONS 



This research was conducted to verify the theoretical 
prediction of shock radius as a function of time. In all 
but four cases, the data collected utilized projectiles im- 
pacting the fluid through a pre-punched hole in the front 
wall of the tank. There were four data points taken where 
the projectile passed through a solid front wall before com- 
ing in contact with the fluid. This was done at the conclu- 
sion of the experiment to see if there was a departure from 
the trends of the previous data and to obtain a qualitative 
estimate of the extent of shock phase damage. The solid 
front wall data were taken using the highest energy level 
tested (12323 in-lb) . 

A. SHOCK RADIUS VERSUS TIME - PRE-PUNCHED FRONT WALL 

Figures 10, 11, and 12 show the experimental data and 
the theoretical prediction of shock radius versus time. 
Seldom were projectile velocities prior to impact constant 
during the testing. Changes in velocity have a direct in- 
fluence on the initial energy. The data were plotted and 
compared, however, assuming constant impact energies. It 
was therefore necessary to determine the shock radius sensi- 
tivity to changes in energy. A plot, using theoretical pre- 
dictions of the shock radius versus impact energy for a 
constant time (figure 13) , indicates that there is an insig- 
nificant change in shock radius for typical experimental 
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Theoretical 
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Figure 10. Shock Radius vs. Time, E_ = 1600 in-lb 



Theoretical 
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Figure 11. Shock Radius vs. Time, E_ = 7493 in-lb 



Solid Front Wall 
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Figure 12. Shock Radius vs. Time, Eq = 12323 in-lb 
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Figure 13. Shock Radius vs. Impact Energy 



variations of impact energy about the reference energy. The 
impact energy calculated from experimental data of the indi- 
vidual tests was within ±200 in-lb of the comparison energy 
level. Corresponding changes in shock radius were of an 
order of magnitude smaller than the accuracy of the radius 
measurements. Therefore, no corrections were applied to the 
measured shock radii for experimental differences in impact 
energy about the comparison energy level. 

The data for an initial impact energy of 1600 in-lb had 
the best agreement with theory. Theory assumes that the 
total projectile energy is deposited at the point of impact. 
The shadowgraphs indicate, however, that none of the projec- 
tiles came to rest at the point of impact. At the higher 
energy levels (7493 in-lb and 12323 in-lb) the data fell be- 
low the theoretical results which implies generation of 
weaker shocks than those predicted by theory. This may be 
explained by considering the projectile as a piston in the 
early stages of shock formation. It would be expected that 
the rate of change of fluid volume displaced by the piston 
determines the momentum exchange between the projectile and 
fluid, hence the shock strength. The shape of the .22 cali- 
ber projectile is noticeably more blunt than the .222 cali- 
ber projectiles (see figure 2). This means that the .22 
caliber projectile would have a greater rate of change of 
fluid volume displaced in the formation stage of the shock 
phase. This would produce conditions which are more closely 
represented by the assumptions of the theory. It is expected 
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that blunt projectiles or warhead fragments produce shocks 
that agree closely with the theory. 

Eventually the shock becomes acoustic because of geomet- 
rical attenuation. The theoretical predictions for the time 
at which the shock is acoustic were 22 micro-seconds at 1.71 
inches for the energy level of 1600 in-lb, 36 micro-seconds 
at 2.81 inches for the energy level of 7493 in-lb and 44 
micro-seconds at 5.40 inches for the energy level of 12323 
in-lb. Theory assumed no counterpressure (Pq<<P) . For this 
reason it is expected that the shock will actually become 
acoustic earlier than theory predicts. The experimental re- 
sults for this study exhibited too much scatter to accurately 
determine when the shock becomes acoustic. It is reasonable 
to assume the duration of the shock phase is accurate in 
order of magnitude. The cavity begins to oscillate at a time 
on the order of 1000 micro-seconds, which implies a distinct 
time separation between the shock and cavity phases. This 
research has shown that the shock phase and cavity phase may 
be studied independently. 

B. SHOCK RADIUS VERSUS TIME - SOLID FRONT WALL 

The four data points taken for the case of a projectile 
passing through a solid front wall and then into the fluid 
are shown in figure 12 as X symbols. Time limitations pre- 
vented the collection of more data. The limited results do 
not indicate a significant departure from the pre-punched 
data. Figures 14, 15 and 16 show typical damage to the .043 
inch thick, 6061-T6 Aluminum front wall specimen. The front 
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Figure 14. Solid Front Wall After Test (6061-T6 Aluminum, 
.043 in. thick, Eq = 12323 in-lb) 




Figure 15. Typical Damage to Solid Front Wall at Point of 
Impact (Front view, 6061-T6 Aluminum, .043 
thick, Eq = 12323 in-lb) 



38 




39 



Figure 16. Damage to Solid Front Wall (Mounted on Tank) 
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wall experienced plastic deformation, became concave in the 
outward direction and buckled. There was some additional 
deformation at a few of the restraining bolt locations be- 
lieved to have been caused by a slightly unsymmetrical load 
distribution. This was a result of projectile impact not 
being at the true center of the front wall. Projectile pene- 
tration caused no visible signs of cracking? although there 
was permanent deformation, the pressure load on the front 
wall was not sufficiently severe to produce cracking. 

Soper {Ref. 7] calculated the percent energy change ver- 
sus impact velocity for penetration of various thickness alu- 
minum plates with no fluid foundation. For the specimen and 
impact velocity used in this study there was just under a 10% 
loss of energy predicted in penetrating the front wall. Using 
the theoretical prediction of shock radius versus impact en- 
ergy in figure 13, there would be an insignificant change in 
the shock radius for such a 10% change of impact energy. 

The shape of the projectiles after penetrating a solid 
front wall was significantly different than the projectiles 
which passed through a pre-punched hole. This was the re- 
sult of the copper jacket on the projectile being deformed 
by both penetration and the stagnation pressure loading on 
the projectile. 

C. PROJECTILE VELOCITY DECAY IN THE FLUID 

The impetus for investigating the projectile velocity de- 
cay in the fluid is the requirement to know the time rate of 
change of energy exchange in order to study the cavity phase 
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[Ref. 1]. It is difficult to predict velocity decay because 
of the many dependent variables of the problem. The analy- 
sis of velocity decay suggested by Lundstrom [Ref. 1] was 
utilized. The equation of motion used was 

D = -mg^ = (4.1) 

dt 2 ^ ^ 

Assuming a constant deformed projectile area, constant drag 
coefficient and constant projectile mass it was possible to 
rearrange equation (4.1) and integrate, which yields 

„ -(p^C^.A/2m)x 

3^ = e ^ ^ (4.2) 

Further rearranging and integration from initial conditions 
of x=0 at t=0 yields the following expression for projectile 
position as a function of time 

X = j2n(BVQt+l) (4.3) 

where 

Pf^D^ 

^ ~ 2m 

The density for the undisturbed fluid, H^O, was assumed con- 
stant. In calculating the theoretical results of equation 
(4.3) the value used for was the projectile velocity im- 
mediately prior to impact with the fluid. It may be more 
suitable to use the velocity of the projectile at the comple- 
tion of shock generation phase as the initial velocity in the 
fluid, Vq. Figure 17 is provided to point out a possible 
discontinuity of the projectile energy and velocity profile 



41 



Projectile Path 



O 

o 

Q> 

Q 







>« 

O 

u 

o 

u 




42 



Figure 17. Qualitative View of Projectile Energy and Velocity Decay in the Fluid 
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during shock formation. The initial location and time of 
this event would still be adequately described by x=0 and 
t=0. A value for the change in projectile energy due to 
shock formation, AE , or the change in projectile velocity 
due to shock formation, AV^, would make it possible to solve 
for Vq. If the experimental data had less scatter than in 
this study it would have been possible to determine the value 
of AE„ required for the theoretical results to be coincident 
with such experimental results. It would then be possible 
to solve for V^. 

Typical projectile shapes after coming to rest in the 
fluid are shown in figure 18 for each energy level tested. 
Figure 4 shows the values for the average final deformed 
area, drag coefficient and mass used to calculate the theo- 
retical curves in figures 19, 20 and 21. Holm [Ref. 8] has 
shown that the projectiles possess dynamic stability and do 
not tumble, which contributes to the deformation experienced. 
Stagnation pressures on the order of 10‘*psi were found to 
exist on the projectile surface. These pressures are large 
enough to cause plastic deformation of the lead .22 caliber 
projectiles and the copper jacketed, .222 caliber projectiles. 

Approximately fifty-one percent (51%) of the .22 caliber 

V 

projectile impact . velocity was lost within 1.5 inches of the 
front wall. The .222 caliber projectiles with an initial 
impact energy of 7493 in- lb only lost 17% of their impact ve- 
locity at 1.5 inches. Projectiles with an energy 12323 in-lb 
lost 25% at the same penetration depth. Percentage differ- 
ences of the energy lost at the same distance from impact 
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Figure 18. Typical Deformation of Projectiles After 
Coming to Rest in the Fluid Filled Tank - 
Pre-punched Hole Results (L. to R. : .22 

caliber, = 1600 in-lb; .222 caliber, 

E = 7493 in-lb; .222 caliber, Eq = 12323 
in-lb) 
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Figure 19 . Projectile Distance Traveled in Fluid vs. Time, Eq = 1600 in-lb 



Theoretical 
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Figure 20. Projectile Distance in Fluid vs. Time/ = 7493 in-lb 
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Figure 21. Projectile Distance in Fluid vs. Time, E = 12323 in-lb 



may be partially attributed to the rate of deformation of 
the projectiles. It is reasonable that the lead .22 caliber 
projectile deforms more rapidly than the .222 caliber jac- 
keted projectile. It appears that the highest energy (12323 
in-lb) projectile has reached a limit at which the projectile 
jacket and lead core cannot sustain the enormous stagnation 
pressure. No longer was the deformation similar to the two 
lower energy projectiles (figure 18). The drag coefficient/ 
deformed projectile area and percentage of projectile veloc- 
ity decay at a fixed distance for the two lower energy levels 
appears consistent. At the highest energy level, values for 
these variables do not appear consistent. When there is ex- 
tensive deformation, as in the case of the highest energy 
level projectile, the assumptions appear to be overly simpli- 
fied. 

The concept of an equivalent flat plate area, f=Cj^A, may 
be more advantageous than considering the drag coefficient 
and deformed projectile area separately. Also included in 
figure 4 is the value for the equivalent flat plate area, f, 
which best matches theory to experiment. This parameter re- 
duces the number of dependent variables by one and describes 
the projectile drag characteristics with one number. 

In reference 7, minimum impact velocities to penetrate 
various thickness plates are given. The experimental veloc- 
ity decay data were extrapolated by theory to predict pro- 
jectile velocities at the rear wall of the tank. These were 
compared with the minimum impact velocity to penetrate the 
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rear wall. After traveling the depth of the tank, 17 inches, 
the .22 caliber projectile (E^ = 1600 in-lb) was predicted 
to have a velocity of 0.394 ft/sec. Experimentally, the .22 
caliber projectiles were found approximately six inches from 
the front wall where equation (4.2) predicts a velocity of 
71 ft/sec. The .222 caliber projectiles with an energy of 
7493 in-lb were predicted to have a velocity at the rear wall 
of 290 ft/sec, which was 32% of the minimum impact velocity 
required for penetration. The experiments with the .222 cal- 
iber projectiles at an energy of 7493 in-lb were the only 
test cases which made a noticeable dent but no penetration 
of the rear wall. The highest energy level projectiles were 
predicted to have a velocity of 105 ft/sec and there were no 
signs of projectile impact with the rear wall. These results 
seem to indicate that the exponential velocity decay is rea- 
sonable where projectile deformation is not extensive. 

D. PRESSURE VERSUS RADIUS AND TIME 

Figures 22, 23, 24 and 25 show theoretical predictions 
of the shock phase pressure pulse and the pressure profiles 
behind the shock at various times. Figure 22 is provided to 
demonstrate the order of magnitude of the shock pressure from 
initial impact until the shock becomes acoustic for an energy 
level of 1600 in-lb. The remainder of the pressure versus 
radius curves were limited to a range of 5 micro-seconds to 
30 micro-seconds. This information was used to estimate 
pressure conditions which will be encountered when further 
experimental studies are undertaken to verify the theory. 
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Pressure - P xlO“® (psi) 



4 - t = Oysec 




E = 1600 in-lb 
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Figure 22. Pressure vs. Radius from Impact - 
Complete Profile, = 1600 in-lb 
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Pressure - P xlO~^ (psi) 




Figure 23. Pressure vs. Radius from Impact, 

E = 1600 in-lb 
o 
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Pressure - P xlO"^ (psi) 




Figure 24. Pressure vs. Radius from Impact, 
Eq = 7493 in-lb 
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Pressure - P (psi) 




Figure 25. Pressure vs. Radius from Impact, 

= 12323 in-lb 
o 
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The pressure profiles behind the shock were not plotted at 
negative pressures although negative pressures behind the 
shock front were predicted in the theoretical results. Since 
water is nearly incompressible, cavitation would occur as the 
pressures became negative. Cavitation has been observed on 
film records of similar hydraulic ram studies done at the 
Naval Weapons Test Center, China Lake and provided for the 
author by Lundstrom. 

Immersion of pressure transducers in the fluid is being 
considered as a means of measuring pressures. Figures 26, 

27 and 28 were plotted to estimate the magnitude of the pres- 
sures experienced at various distances from the front wall 
and the duration of the pressure pulse. The plots indicate 
the need for a transducer with a very fast response time. 

It is expected that actual pressures will be somewhat 
lower than predicted, with the possible exception of blunt 
nose projectile data. There are two reasons for predicting 
lower pressure profiles. As explained earlier, the shock 
radius versus time data implied a weaker shock than the theo- 
retical predictions. This was the result of the assumption 
in the theory that total impact energy was deposited in the 
fluid during shock formation. A weaker shock would result 
in correspondingly weaker shock pressures throughout the 
shock phase. The second reason pertains to the pressure pro- 
file as the shock Mach number approaches one. At lower shock 
Mach numbers, the assumptions of undisturbed fluid pressure 
in front of the shock being much less than at the shock 
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Figure 26 . Pressure vs. Time for Three Radii, Eq = 1600 in-lb 
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Figure 27. Pressure vs. Time for Three Radii, E = 7493 in-lb 
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Figure 28. Pressure vs. Time for Three Radii, E = 12323 in-lb 



CPq<<P} will become invalid. This assumption should cause 
the theoretical predictions to be higher than experimental 
results. 
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V. RECOMMENDATIONS 



1. The scatter of the experimental data measuring shock 
radius versus time could be reduced by placing the fourth 
chronograph screen closer to the tank. New counters with 
high frequency, temperature compensated oscillators to 
provide a better time base are also recommended. 

2. The ballistic range should be modified to provide the 
capability of firing higher energy projectiles. Shock 
radius versus time and projectile distance traveled in 
the fluid versus time should be studied at these higher 
energy levels. 

3. Fluid pressure at and behind the shock should be com- 
pared with theory at various energy levels. 

4. This study used lead or copper jacketed, lead core 
projectiles and no tumbling was observed. Further stud- 
ies using steel jacketed projectiles should be completed 
to test for deviations from the data of this study. This 
information would provide the fuel cell designer with 
data typical of ballistic threats encountered in combat. 

5. Strain gauges should be placed on the front wall spec- 
imens to determine the stresses incurred during the shock 
phase. Strain data could then be compared with SATANS 
computer code . 

6. Higher energy projectiles which produce cracking or 
catastrophic failure of the front wall specimen should 
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be tested. These data would provide the designer with 
the critical energy level at which shock phase damage 
may result in an aircraft loss. 
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COMPUTER PROGRAM 



* This PROGRAM CALCULATES THE FOLLOWING 

* PART .1 

* 1. SHOCK RADIUS VERSES TIME 

* 2. SHOCK MACH NUMBER VERSES TIME 

* PART II " 

« 1. PRESSURE VERSE TIME AND RADIUS 

* FROM IMPACT BEHINC THE SHOCK 

* FOR A GIVEN FLUID AND IMPACT ENERGY, EG. THIS PROGRAM 
« ASSUMES: SHOCK RADIUS IS PROPORTIONAL TO TIME TO THE 

* C.8 POWER UNTIL ACOUSTIC, A POWER LAW DENSITY PROFILE 

* BEHIND THE SHOCK, A STRONG SHOCK, ADIABATIC ACROSS 

* THE SHOCK, THE UNDISTURBED FLUID PRESSURE IS ESSEN- 

* TIALLY ZERO, SPECIFIC INTERNAL ENERGY IN UNDISTURBED 
« FLUID IS ESSENTIALLY ZERO. 



❖ ^ 35; :{s ;{£ :{s :{c :-c 4c 3{c :jc 3}C :{s ajc :jc :-c :-c ^ :jc :-c :{: :{s 



SYMBOL DESIGNATION 

EO= IMIPACT ENERGY OF PROJECTILE (IN-LEJ 
CO= SPEED OF SOUND FOR FLUID (FPS) 

RHOO= DENSITY IN UNDISTURBED FLUID (LE/FT3) 

DC= INCREMENT OF NON-DIMENSIONAL RADIUS, DELTA 
ZHETA 

DELT= INCREMENT OF TIME (MICRO-SEC) 

ET= ENERGY/VOL. = 2-EO/ (R*«3 ) 

Fl= NON-DIMENSIONAL PRESSURE AT THE SHOCK VS. 
SHOCK MACH NO. 

DODM= DERIVATIVE OF NON-DIMENSIONAL PARTICLE 
VELOCITY W.R.T. SHOCK MACH NO. VS SHOCK 
MACH NO. 

VS= SHOCK MACH NO. 

X= PERCENT EG LOST AT ENTRY WALL 



DIMENSION ET(20) , PSI 1(20) ,Fi(20) ,VS(20) , DCDM(20) 
1,RXMS(20) ,DRS(20), ARMS(20) ,DT(20),T(20) ,RE(20) , 
1RV(20) , RDCD(20) ,XMS(20) ,RS(20) , RPS I (20 ) , R F ( 20 ) 
DIMENSION TT(300) ,XM( 300) , R ( 30C ) , F2 ( 300) , PS 1 2 (300 ) , 
iDCDM'2(300) ,Q(300) ,C(100) ,F(2 ,20 ) ,P(2 ,20) 

DIMENSION CCNRS(30C) ,SHGCV(300) ,TCFGR(300) 

DIMENSION F0R(20) 

DC 500 J=i,20 
FOR(J)=C.G 
5C0 CONTINUE 
N = 18 

READ (5,100) EO,CO,RHOO,X,DC,DELT 
READ (5,101) (ET( I ) , I = 1,N ) 

READ ( 5, 101) ( PSI 1 ( I ) ,I = 1,N) 

READ ( 5 , 101) ( Fl( I ) , 1 = 1, N ) 

READ (5,101) (VS( I) , 1=1, N ) 

READ (5,101)(00DM( I) ,I = 1,'\) 

WRITE (6,200 ) EO,CO,RHCO,X,DC,DELT 
DC 15 1=1, N 

15 WRITE (6,201) ET ( I ) , P S 1 1 ( I ) , F 1 ( I ) , VS( I ) , DOOM ( I ) 

PART I 
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PERCENTAGt OF INITIAL ENERGY DEPOSITED ANC INITIAL 
SHOCK RADIUS 

E0=X<'E0 

R0=(E0<'32..2=S'144./ (RH00*C0*C0*2.*3.1416) ) *=t0. 33333 
INITIAL VALUE OF ENERGY/VOL. 

ET0=2.=«=EC/(R0*<'3.0) 

INVERT INPUT ARRAYS FOR USE WITH INTERPOLATION 
ROUTINE, PIF2 

DC 12 1=1, N 

J=N-I+1 

PEU)=ET(J) 

RPSKI )=PSI1(J) 

RFCIJ=F1(J) 

RV( I) =VS( J) 

12 RCCDd )=00DM(J) 

INTERPOLATE FOR INITIAL VALUES AT SHOCK 

VSO=PIF2(ETO,RE,N,RV) 

F10=PIF2(VS0,RV,N,RF) 

PSI10 = P IF2(VS0,RV,N,RPSI) 

DCDK0=PIF2 (VSO,RV,N,RDCC) 

COMPUTE INITIAL TIME FOR DISTANCE RO 

TC = R0/(CC*VS0=5'12.0) 

ASSIGN SUBSEQUENT VALUES OF SHOCK MACH NC . AT WHICH 
SHOCK RADIUS WILL BE COMPUTED 

DC 10 I=1,N 

IF (VSO-VSm) 1,2,2 

1 GO TO 10 

2 M=N-I+1 
XMS (1 )=VS0 
GO TO 11 

1C CONTINUE 

COMPUTE SUBSEQUENT VALUES OF SHOCK RADIUS 

11 RS(i)=R0 

DO 20 J=2,M 

XM.S(J)=VS(I) 

RS( J) = (2.=S'E0 /(ET(D) )*«0.3333 
1=1 + 1 

20 CONTINUE 
K=M-1 

COMPUTE DELTA TIME CORRESPONDING TO SHOCK MACH NO. 
AND SHOCK RADIUS 

DC 30 1=1, K 

DRS(I)= RS(I+1)-RS (I) 

ARMS! I )=2-/(XMS( I )+XMS( 1+1} ) 

30 OT ( I) =ORS ( I) =^ARMS( I) / (00*12.0) 

T( 1) = T0 
DC 40 J=2,M 

40. T(J)= T( J-l)+OT( J-1) 

WRITE (6,102) 

WRITE (6,103) EO 
WRITE (6,110) 

PRINT SHOCK MACH NC. AND SHOCK RADIUS VS TIME 
DC 50 1=1, M 

50 WRITE (6,111) T( I ) ,XMS( I ) ,RS ( I ) 
lOOFORMAT (6E10.0) 
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ICl FGRNAT (0E1O.O) 

102 FORMAT (1H,//,39X, • ENERGY 15X , 'FLUID* ) 

103 FORMAT (29X,1E20.7) 

lie FORMAT (1H,//,39X • TI ME* , 17X MACF ' »17X RADIUS • ) 

111 FORMAT (29X,3E20.7J 
2CC FORMAT (//, 6E20.7) 

2C1 FORMAT (15X,5E20.7) 



PART II 



INITIALIZE VARIABLES AND ASSIGN VALUES USED WITH 
NON-DIMENSIONAL PRESSURE DISTRIBUTION FUNCTION 

BEB=T(M)*l.E0e+1.05 
L = BBB 
TT(1)=0.0 
XM( 1)=XMS(1) 

R(1)=R0 

F2{1)=F10 

PSI2(1 )=PSI10 

DCDM2(1 )=DCOMO 

Q(1)=3.=:=(PSI10-1.0) 

e=-0.25 

1=1. 0/DC +1.05 
C(1)=0.C 
WRITE (6,301) 

WRITE (6,302) TT ( 1 ) ,XM (1 ) , R (1 ) 

WRITE (6,303) 

DC 60 J=2,N 

6C C(.J)=C( J-D+OC 

CCM.FUTE NON-DIMENSICNAL AND CIMENSICNAL PRESSURE FOR 
INITIAL SHOCK RADIUS. 

DC 70 J=1,I 

F(1,J) = F2(1)+PSI2(1)«(C(J)*«(Q(1)+2.0)-1.C)=^(F2(1)/ 
lPSI2(l)-e^(F2(l)+XM(l)=i=CG0M2(i) ) ) / ( Q ( 1 ) + 2 .0 ) 

P(1,JJ=F( 1 tJ)-RHCO^CO*CO*XM(l)=PXM(l)/(32 .2=^144. ) 

PUNCH 91,R(1),C(J) ,P(1,J) 

70 CONTINUE 
K=1 
JJ=1 
SUM=0.0 

CALL FORCE(FOR,C,P,R,I,JJ,K) 

DC 501 J=1,I 
SLM=SUM+FGR( J) 

WRITE (6, 304) C( J ) , P ( 1 , J ) , F ( 1, J ) , FCR( J ) 

5C1 CONTINUE 

TCFCR( 1)=SUM 

ITERATION TO SOLVE PRESSURE, SHOCK RADIUS, SHOCK MACH 
NO. VS TIME 

DC 80 K=2,L 

TT(K)=TT(K-1)+DELT 

IF (TT(K)-TO) 21,21,22 

21 XM(K)=XM(1) 

R(K)=R( 1) 

F2(K)=F2(1) 

PSI2(K)=PSI2(1) 

DCDM2(K)=D0DM2( 1) 

G(K)= Q(l) 

DC 55 J=1,I 

55 P(2,J)=F(1,J) 

PUNCH 92,R(K) 

GC TO 80 

22 X=TT(K) 

INTERPOLATION OF INPUT DATA AND PART I DATA FOR 
VARIABLES IN PRESSURE FUNCTION 
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XK(K)=PIF2(X,T,M,XNS) 

R(K)=PIF2(X,T,M,RS) 

Y=XM(K) 

F2(K)=PIF2(Y,RV,N,PF) 

FSI2(K)=PIF2(Y,RV,N,RPSI) 

DC0M2(K)=PIF2 (Y.RV,K,ROCD) 

Q(K)=3.=fr(PSI2(K}-i.G) 

WRITE (6,301) 

WRITE (6,302) TT ( K ) , XP ( K ) , R ( K ) 

WRITE (6,303) 

DC 90 J=1,I 

F(2,J)=F2(K)+PSI2(K)*(C(J)=i'«(Q(K)+2.0)-1.0)^(F2(K)/ 
1PSI2(K)-8«(F2(K)+XM(K )*D0DM2(K) ))/(Q(K)+2.0) 
P(2,J)=F(2,J)*RH0G=«=C0*C0=i=XM(K)*XP ( K ) / ( 32 . 2* 144 . ) 

PUNCH 91,R(K) ,C(J) ,P(2, J) 

90 CGNTINUE 
JJ = 2 
SUN=0.0 

CALL FCRCE(FCR,C,P,R,I,JJ,K) 

DC 502 J=1,I 
SLM=SUM+F0R(J) 

WRITE (6 ,3X)4) C( J) , P ( 2 , J ) , F ( 2 , J ) , FCR ( J ) 

502 CONTINUE 

TCFCR(K)=SU.M 
8C CONTINUE 

WFITE(6,305) 

CALL CCNVER(R,XM,TCFOR,TT,L) 

91 FORMAT (3E15.5) 

92 FORMAT (1E15.5) 

2C1 FORMAT ( IH ,// ,36X, ' TI ME ', 17X, ' MACF ', 17X ,' RADIUS ' ) 

302 FORMAT (29X,3E20.7) 

303 FORMAT (1H,//,39X, • P/RS' ,17X , 'P ' ,17X, • F' , 19X, 'FORCE' ) 

304 FGPMAT(29X,4E20.7) 

305 FORMAT (lh,//,36X, ' T I ME ' , 17X , ' V ELC ' , 17X, ' R ACLI S ' , 1 7X , ' T 
IGFGRCE' ) 

ENC 

FUNCTION PIF2 ( X , XL 1ST , N , FL 1ST ) 

C 

C SECOND ORDER INTERPOLATION 

C 

DIMENSION XLIST (ICO), FLIST (100) 

BLIF (P,Q,R,S,T) = ( (Q-P)*(S-T)/(R-Q)+S) 

IF (X-XLIST(N)) 2,1,1 

1 I = N-1 
GO TO 5 

2 IF(X-XLISTd) ) 4,4,6 

4 1 = 1 

5 K = 1 

GO TO 30 

6 K » 2 

7 DO 8 I = 1,N 

IF (X-XLIST(D) 9,9,8 

8 CONTINUE 
I=N 

9 1 = 1-1 

30 BLIFl = BLIF(X,XLIST(I),XLIST( I+1),FLIST(I),FLIST(I+1) 
1) 

10 IF (K-1) 11,11,12 

11 PIF2 = BLIFl 
RETURN 

12 IF((I+2)-N) 13,13,16 

13 IF ((I-i)-l) 15,14,14 

14‘ IF(A6S( XL IST( I-1)-X)-ABS(XLIST( I + 2)-X) ) 16,15,15 

15 L = 1+2 
GO TO 17 

16 L = I-l 

17 BLIF2 = BLIF (X , XL 1ST ( I ) , XL 1ST ( L ) , FL 1ST ( I ) , FL I ST ( L ) ) 
PIF2 = 6LIF (X,XLIST( I+l) ,XLIST( L) , BLIFl, 6LIF2) 

18 RETURN 
ENC 
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SLEROUTINE FORCE ( FOR ,C , P ,R , I , JJ ,K) 



* 

* 

« 

* 

« 

« 

« 

* 

* 

♦ 

=* 

« 

* 

« 

* 

* 

4 < 



* 

* 

This SUBROUTINE COMPUTES THE FORCE ON THE ENTRY WALL * 
OLE TO SHOCK PHASE FOR A FRONT WALL WITH A « 

PRE-PUNCHEC HOLE, I. e. PROJECTILE COES NOT IMPINGE CN * 
FRONT WALL UPON ENTRY. COMPUTES FORCE TC EE ZERO FCR=> 
PRESSURE FIELD OVER THE PRE-PUNCHEC HOLE CR WHERE THEWS' 
PRESSURE GOES NEGATIVE. THE PROGRAM ACCOMPLISHES THIS=!= 

❖ 
❖ 

* 
* 
* 



BY A STAIGHT LINE APPROXIMATION EETWEEN PRESSURE 
DATA POINTS. 

RHRS= NON-DIMENSIONAL RADIUS OF PRE-PUNCHEC HOLE 
S= SLOPE OF STRAIGHT LINE BETWEEN PRESSURE DATA 
POINTS 

RX= THE INTERCEPT OF R/RS , NON-C IMENS lONAL RADIUS, 

WHERE PRESSURE IS ZERO 
FCR(N)= FORCE(LB) FOR AN ANNULUS WHOSE OUTSIDE RADIUS* 
IS C(N) AND INSIDE RADIUS IS C(N-l) * 

* 






DIMENSION F0R(20) ,C(100) ,P( 2,20) ,R(300) 
RHRS=.375/R(K) 

IF(RHRS-i. ) 100,25,25 
IOC IF(JJ-l) J.,1,2 

1 DC 10 N=2, I 

S=(P( 1,N)-P( 1,N-1) )*4. 

IF(P(1,N-D) 3,4,4 

3 IF(P(1,N)) 5,5,6 

5 FCfi(N)=C.D 
GO TC 10 

6 RX={S*C(N-1)-P(1,N-1) )/S 
IF(RHRS .GT.C(N) ) GO TO 43 
IF (RX. GT.RHRS) GO TC 60 
PH=S*(RHRS-C(N-1) )+P( 1,N-1) 

FCR(N) = (.5*(P(i,N)-PH )+PH)*3.1416*(C(N)*C (N)-RHRS 
1*RHRS)*R(K)*R(K) 

GO TO 1C 

60 FCR(N) = P(1 ,N)*.5*3 .1416*(C (N)*C(N )-RX*RX ) *R (K ) *R ( K ) 
GO TO 10 

4 IF(RHRS-C(N-1) ) 40,40,41 

4C FOR(N) = (.5*(P(1,N)-P( 1 , N-1 ) ) + P ( 1 , N-1 ) ) *3 . 14 16* ( C ( N ) 
1*C(N)-C(N-1)*C(N-1) )*R(K)*R(K) 

GO TO 10 

^1 IF(RHRS-C(N) ) 42,43,43 

42 PH=S*(RHRS-C(N-1) )+P( 1,N-1) 

FCR(N) = (.5*(P(l,N)-Fh)+PH)*3.14ie*(C(N)*C(N)-RHRS 
1*RHRS)*R(K )*R(K) 

GO TO 10 

43 FOR(N)=C.O 
1C CONTINUE 

RETURN 

2 DC 20 N =2 , I 
S=(P(2,N)-P(2,N-1) )*4. 

IF(P(2,N-D) 13,14,14 

13 IF(P(2,N)) 15,15,16 

15 FCR(N)=0.0 
GO TO 20 

16 RX=(S*C(N-1)-P(2,N-1) )/S 
IF(RHRS.GT.C(N) ) GO TO 143 
IF(RX. GT.RHRS) GO TO 160 
PH=S*(RHBS-C(N-i ) ) +P( 2,N-1 ) 

FOR (N)=(.5*(P (2, N)-PF)+FH) *3.1416* (C(N)*C(N)-PHRS 
1*RHRS) *R(K )*RtK) 

GO TO 20 

16C FCR(N) = P(2,N)*.5*3.1416*(C(N)*C(I\)-RX*RX)*R (K)*R(K) 
GO TO 20 

14 IF(RHRS-C(N-1) ) 140,140,141 
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140 FCR(N)=( .5«(P(2tNi-P(2,N-i) )+P( 2,N-1) )*3.1416*(C(N) 
l#C(N)-C(N-i)=frC(N-l) )-R(K)*R(K) 

GC TO 2C 

141 IF(RHRS-C(N) ) 142,143tl43 

142 PH=S*(RhRS-C(N-l) )+P(2,N-i) 

FGR(N) = ( .5'-^(P(2,N)-PH)+PH>*3.141£*(C(N)*C{Ni-RHRS 
1#PHRSJ^P(K)*R(K) 

GO TO 20 

143 FCR(N)=0.0 
2C CCNTINUe 
25 RETURN 

ERC 

SLEROUTINE CCNVER ( R , XM , TOFOR ,TT , L ) 

THIS SUBROUTINE CONVERTS VELOCITY AND RACIUS 
TO KM/SEC AND CM RESPECTIVELY 

CIMENSICN Rt300J ,XM(300),TOFOR(300),CONRS(3C0) ,SHOCV( 
13CC) ,TT(30C) 

DC 510 11=1, L 
CCNRSnn = R(II)=«'2.54 
ShCCV( I I ) = XM( II)={'4.9/3.281 

WRITE (6, 520) TT( II ) , SHuCV ( 1 1 ) , CONRS ( 1 1 ) ,TCFCR(II) 

52C F0RMAT(29X,4E20.7) 

510 CCNTINUE 
RETURN 
END 
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